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surface.  Through  such  basic  research,  knowledge  of 
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contribute  to  advancing  the  state-of-the-art  in 
weather  forecasting. 
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SUMMARY 


The  objective  of  this  study  is  to  develop  a 
theoretical  model  for  the  structure  of  turbulence  in 
the  atmosphere  and  to  solve  the  equations  for  the 
distribution  of  wind  and  turbulence  in  the  planetary 
boundary  layer.  Starting  with  the  basic  equation  of 
motion  for  an  incompressible  fluid,  it  is  modified  to 
incorporate  the  mixing-length  hypothesis  of  Prandtl  to 
relate  the  turbulent  stresses  to  the  mean  flow 
characteristics.  It  is  assumed  the  atmosphere  is 
adiabatic,  barotropic,  and  in  a  steady  state.  These 
assumptions  are  not  all  essential  to  the  solution , but 
do  simplify  the  discussion.  Based  on  the  assumptions, 
a  relation  for  the  mixing-length  distribution  within 
the  boundary  layer  is  developed.  Using  this  relation¬ 
ship  in  the  equation  of  motion  led  to  a  set  of  second 
order,  nonlinear  differential  equations,  which  were 
solved  on  a  digital  computer.  Universal  profiles  of 
the  wind,  stress,  and  eddy  viscosity  were  fixed  by 
invoking  the  important  notion  of  similarity;  that  is, 
it  Is  assumed  the  scale  of  turbulence  is  uniquely 
related  to  the  gross  dimensions  of  the  boundary  layer. 
The  requisite  universal  constant  is  evaluated  from 
experimental  data.  Possible  applications  of  the  model 
to  practical  problems  are  outlined. 
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NUMERICAL  SOLUTION  OF  THE  VISTRI BUTION 
OF  til INV  AMP  TURBULENCE 
IN  THE  PLANETARY  BOUNDARY  LAYER 

PA  Taife  I-A -0-11001-B-021-08 


I.  INTRODUCTION 


Meteorology  has  long  been  divided,  into  the 
subordinate  fields  of  macrometeorology  and  micrometeorology. 
This  division  has  not  been  entirely  arbitrary,  for  the 
interests  and  successes  in  these  two  areas  have  been 
rather  widely  separated.  If  we  were  to  select  the  most 
characteristic  elements  from  each  of  these  divisions, 
the  most  logical  are  the  geostrophic  wind  to  characterize 
macrometeorology  and  the  logarithmic  wind  profile  for 
micrometeorology.  doth  of  these  relationships  are 
rather  successful  formulas.  The  geostrophic  wind  equation 
lias  been  used  to  evaluate  the  pressure  distribution  from 
the  horizontal  field  of  wind  in  the  free  atmosphere  , 
whereas  the  logarithmic  wind  profile  equation  has  been 
used  to  relate  the  vertical  distribution  of  wind  speed 
immediately  above  the  earth's  surface  to  the  momentum 
transport  and  the  roughness  characteristics  of  the 
surface.  On  the  face  of  it,  there  does  not  appear  to  be 
any  direct  connection  between  these  two  relationships 
for  they  deal  with  different  regions  of  the  atmosphere. 
However,  it  is  the  intent  of  this  report  to  present  a 
hypothesis  which  unifies  these  two  concepts  in  such  a 
manner  as  to  predict  the  vertical  distribution  of  wind 
and  turbulence  characteristics  within  a  barotropic, 
adiabatic,  steady-state  planetary  boundary  layer. 
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II.  VE  l/E  LOPMENT  OF  THE  HOVEL 


The  basic  principle  which  is  involved  in  this  study 
-  is  the  equation  of  motion: 


dV  3V  -+  ->  „  1  * 

dt  ~  Tt  +  voVv  --2QXV  +  g  +  -  VOX 

where  v  is  the  vector  wind  velocity  relative  to  a  point 
on  the  earth’s  surface,  n  is  the  angular  velocity  of 
the  earth’s  rotation,  g  lo  the  acceleration  of  gravity, 
and  x  is  the  stress  dyadic.  Equation  (1)  is  notoriously 
difficult  to  solve  without  considerable  simplification. 

For  this  study  a  number  of  simplifying  assumptions 
are  invoked.  We  assume  that  the  wind  is  steady  and 
uniform  in  the  horizontal  and  that  quasi-hydrostatic 
equilibrium  exists  along  the  vertical  axis.  As  a  result, 
equation  (1)  can  be  written  as 


l 

-fu  =  -  3p 
p  3y 


1  Hi 

p 


fv  =  A  iE  _  A  3tx 

p  ax  ”  p  9Z 


(2) 


where  v  =  ui  +  vj  ,  f  =  2Q i-cne  where  0  is  the  latitude, 
p  is  the  pressure,  and  tx  and  iy  are  the  stresses  in 
the  i  and  j  directions,  respectively.  At  this  point 
it  is  possible  to  introduce  the  geostrophic  wind  (G) 


which  is  defined  as 


G  —  Gv  i  +  G, 


Substituting  in  equation  (2)  we  obtain: 

f(Gx-u)  +1  Hz  _  n 
X  o  3z  “  0 


-  (12  t  - 
f  \3y  1 


-f ( G  -v )  +  A 

y  P  3Z 


o 


(3) 


rJ  a, 


With  fev;  exceptions,  the.  atmospheric  boundary 
layer  is  fully  turbulent  and  thus  for  practical  purposes 
the  stresses  of  equation  (3)  are  wholly  the  Reynolds 
stresses.  The  Reynolds  stresses  are  essentially 
statistical  quantities  that  express  the  covariances  of 
the  turbulent  fluctuations.  The  solution  of  equation  (3) 
thus  requires  a  knowledge  of  the  statistical  functions 
describing  everywhere  the  turbulent  motions.  Although 
there  have  been  in  recent  years  significant  advances  in 
understanding  the  nature  of  turbulence,  our  knowledge 
still  remains  inadequate  to  derive  transport  relation¬ 
ships  from  the  primitive  equations.  Consequently, 

Imost  every  attempt  to  solve  turbulent  transport 
roblems  must  lean  heavily  on  semi-empirical  or  pheno¬ 
menological  theories.  The  "mixing-length'’  hypothesis 
of  Prandtl  has  played  the  leading  role  in  the  area  for 
three  decades.  This  concept  continues  in  one  way  or 
another  to  be  used  in  turbulent  transport  models  simply 
because  it  gives  useful  results.  For  this  very  reason, 
we  shall  utilize  the  mixing-length  hypothesis  to  solve 
equation  (3)  knowing  full  well  that  the  results  will 
not  be  precise,  but  yet  believing  that  they  will  be 
useful. 


For  our  purposes,  the  mixing-length  relationship 
cat;  be  expressed  by 

t/o  =  U*  =  (i  3U  )2 

“  (*Q 

wnere  U#  =  (t/o)i/2  is  the  friction  velocity  and  l  is  the 
mixing  length.  This  expression  is  more  commonly  seen 
in  the  form: 


t  /  o  -  Ugb  9_U  _  r  5U 
3z  3z 


where  K  is  the  eddy  viscosity.  There  have  been 
conflicting  views  as  to  what  factors  control  the  magnitude 
of  the  mixing  length.  .'on  Karman  took  the  view  that  the 
magnitude  cf  the  mixing  length  is  determined  by  local  flow 
conditions.  The  von  Karman  hypothesis  has  certain 
undesirable  featares  and  in  general,  has  not  been  accepted, 
(dee  Lettau1).  On  tne  other  hand,  the  notion  of  Prandtl 
has  generally  been  that  the  magnitude  of • the  mixing 
length  depends  upon  gross  features  of  the  flow.  In  the 
case  of  free  turbulence,  the  mixing  length  would  be 
assumed  constant  in  a  cross  section  of  the  mixing  zone 


i 


and  its  magnitude  proportional  to  the  width  of  the 
mixing  zone.  On  the  other  hand,  in  a  wall  turbulence 
regime  it  is  assumed  that  the  magnitude  of  the  mixing 
length  is  proportional  to  the  distance  from  the  wall. 
For  our  purposes,  this  can  be  expressed  by 

l  -  kz  (6) 

A  A 

where  k  is  von  Harman* s  constant.  Using  eguations  (4) 
and  (6),  the  logarithmic  wind  profile  for  the  surface 
boundary  layer  can  be  derived  by  integration,  that  is 

U  =  U*ln(z/z  ) 

£  °  (7) 

« 

where  zq  is  the  roughness  length. 

Neglecting  the  variation  of  density  with  height, 
we  can  incorporate  equation  (5)  into  equation  (3)  with 
the  result: 

f (G  -u)  +  U,  L  av 

3 z  3z  "  0 

-f(Gy-v)  +  l_  U,  £  3u  ,  0 


In  actuality,  equation  (7)  is  a  contradiction  to 
equation  (8).  The  logarithmic  wind  profile  is  based 
on  the  notion  that  the  shearing  stress  U#  l  du  is 

37 

constant  with  height, whereas  equations  (8)  clearly 
reveal  that  the  vertical  gradient  of  sheading  stress 
is  actually  greatest  near  the  boundary.  This  paradox 
can  be  reconciled  by  recognizing  that  the  percentage 
rate  of  change  of  stress  near  the  boundary  is  so  small 
that  for  practical  purposes  the  velocity  profile  is 
essentially  logarithmic.  At  the  other  extreme  of  great 
height,  the  stress  will  vanish  and  by  equation  (8) 

*  A  more  general  discussion  of  the  roughness  length 
can  be  found  in  ERDAA-MET-7-63  report,  **A  Model  for 
Wind  Flov  in  an  Idealized  Vegetative  Canopy.” 
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the  wind  becomes  geostrophic,  thus  the  solution  of 
equations  (8)  should  possess  the  requisite  features 
desired  in  our  objective. 

The  most  notable  early  attempt  to  solve  equations 
(8)  with  variable  eddy  viscosity  was  that  of  Rossby 
and  Montgomery2.  They  proposed  a  two-layer  model. 

In  a  relatively  shallow  layer  near  the  ground  they 
assumed  the  mixing  length  was  linearly  increasing  as 
indicated  by  equation  (6).  The  eddy  viscosity  was 
assumed  to  be  a  maximum  at  the  top  of  this  layer. 

Above  this  was  a  relatively  thick  transition  layer  in 
which  the  eddy  viscosity  decreases  as  a  quadratic 
function  of  height.  The  model  of  Rossby  and  Montgomery 
was  a  definite  improvement  over  previous  models  which 
assumed  a  constant  eddy  viscosity.  The  major  criticism 
of  the  Rossby-Montgomery  model  is  that  the  manner  in 
which  the  eddy  viscosity  is  specified  is  quite  arbitrary. 
It  is  the  intent  of  this  study  to  obtain  unified 
solutions  of  equations  (8)  on  the  basis  of  an  assumed 
mixing- length  distribution.  It  is  proposed  that  on 
the  basis  of  similarity  the  magnitude  of  the  mixing 
length  is  determined  by  gross  features  of  the  flow. 
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NUMERICAL  SOLUTION  OF  PI  FFERENTI AL  EQUATIONS 


According  to  equation  (5)  the  eddy  viscosity  is 
equal  to  the  product  of  the  friction  velocity  and  the 
mixing  length,  and  according  to  equation  (4)  the 
friction  velocity  depends  on  the  mixing  length  and  the 
wind  shear.  In  the  surface  boundary  layer  the  friction 
velocity  behaves  as  if  it  were  nearly  independent  of 
the  height,  whereas  at  great  height  the  wind  approaches 
the  geostrophic  wind  and  the  friction  velocity  becomes 
negligible.  On  the  other  hand,  the  mixing  length, 
according  to  equation  (6)  should  be  valid  ad  infinitum.. 
The  outer  portion  of  the  boundary  layer  behaves  more 
like  a  region  of  free  turbulence  than  of  wall  turbulence, 
and  consequently,  one  would  be  led  to  propose  that  the 
mixing  length  approaches  a  limiting  value  In  the  layers 
remote  from  the  boundary. 


In  this  study,  it  is  proposed  that  when  the 
planetary  boundary  layer  is  in  an  adiabatic  steady 
state,  the  mixing  length  increases  with  height  at  the 
rate  given  by  von  Karman's  constant,  but  that  there  is 
a  direct  linear  feedback  which  prevents  unbounded 
growth.  This  notion  is  expressed  by  the  following 
differential  equation: 


d l 
dz 


■  k  -  l/L 


(9) 


where  L  Is  a  constant. 

The  solution  of  equation  (9.)  statisfying  the  boundary 
conditions  is: 


l  -  l  (l-e“z/L) 
o 


(10) 


where  £  /L  »  k.  This  equation  has  the  characteristics 
that  l  essentially  increases  linearly  with  height  near 
the  boundary, but  approaches  the  limiting  value  l0  at 
great  height.  Figure  1  compares  equation  (1)  with  the 
mixing-length  distribution  as  computed  from  Lettau's3 
reanalysis  of  the  Leipzig  wind  profile.  Equation  (10) 
fits  the  data  reasonably  well  except  in  the  upper 
portion  of  the  profile.  However,  in  this  region  the 
evaluation  of  the  mixing  length  Is  much  less  reliable 
because  of  the  small  values  of  the  velocity  derivative. 
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For  this  study  equation  (10)  has  been  utilized  as 
a  general  expression  for  tne  mixing  length.  Having 
this  express  i-m  it  is  possible  to  solve  for  the 
velocity  profile  in  the  planetary  boundary  layer.  For 
computational  purposes,  it  is  necessary  to  express  the 
equation  in  nondimensional  form.  For  this  purpose  we 
define  the  following  variables: 

U  =  (Gx-u)/G,  V  =  (Gy-v)/G,  s—  tx/PG2,  t=  ry/pG2, 

T  =  U_/G  ■  (s2  +  t2)1/*4  R  =  G/fz.  x  =  z/z_ 

*  o  o  o 

wnere  R  is  termed  the  Rossby  number.  Using  these 
expressions,  equations  (3)  can  be  expressed  as 


V  =  12.  U 

°  3x  , 


(11) 


and  equation  (5)  can  be  expressed  as 

s  =  -  T<VZo)r  If  t  -  -  ?U0/z0)y  H 

’  (12) 
where  .y  =  1  -  exp(-z/L).  If  we  differentiate 
equations  (11)  and  substitute  the  results  in  equations 
(12),  we  obtain 

s  -  h  Ty  ill  t  =  -  AnTy  a2s 

T5T7  (13) 


where  AQ  =  RQ  tQ/zQ.  Equations  (13)  are  the  basic 
equations  to  be  solved  for  the  boundary  layer  model. 

The  boundary  condition  requires  that  s  *  t  *  o  at 
great  height  and  we  chose  axes  of  orientation  such  that 
at  the  lower  boundary  t  =  o  and  s  -  s0.  For  each 
value  of  the  parameter  A0,  there  must  be  a  corresponding 
value  of  sQ,  However,  the  actual  relationship  of  s0 
to  A0  cannot  be  determined  until  the  solutions  are 
known.  To  obviate  this  paradox  it  is  necessary  to 
rewrite  equations  (13)  in  a  form  containing  only  a 
single  parameter.  This  was  done  by  defining  new 
variables  as  follows: 
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Utilizing  these  terms,  equations  (13)  become 


Q  =  10 1 2  Py 


32H 


ax' 


R  =  -  10 


1  2 


a  2q 

p’' 


OH) 


with  the  following  boundary  conditions 

x  -  1  Q  =  Q  R  =  0 

o 

X  +  ”  Q  =  R  a  p  =  0 

Prom  equation  (I1!)  and  the  boundary  conditions,  it 
can  be  seen  that  a  particular  solution  is  completely 
specified  by  the  quantity  Q0  provided  y  is  solely  a 
function  of  proximity  to  the  boundary. 

Considerable  difficulty  was  encountered  in  solving 
equations  (14).  It  was  necessary  to  assume  that  Q  and 
R  became  negligible  at  a  finite  but  large  height  H 
rather  than  at  infinity.  Numerous  mathematical  methods 
were  tried  to  solve  the  equations,  but  the  solutions 
were  finally  achieved  more  by  utilizing  physical 
intuition  than  mathematical  prowess.  Briefly,  the 
method  of  solution  consisted  of  transforming  the 
independent  variable  to  y  =  2.0  x  and  then  establish¬ 
ing  a  grid  of  100  equally  spaced  points  in  the  y  domain 
between  x  =  1  and  x  =  H,  The  derivatives  at  the  grid 
points  were  approximated  from  second-order  Laprangian 
polynomials  utilizing  the  values  of  the  variables  at 
adjacent  grid  points.  On  the  basis  of  the  ordinary 
differential  equations,  a  form  of  parabolic  partial 
differential  equation  was  assumed  and  the  final  solution 
was  derived  through  integration  until  the  solution  of 
the  partial  differential  equations  decayed  to  the  steady 
state  ordinary  solution. 

Choosing  z  /£0  a  3.125  x  lo”5,  numerous  solutions 
of  equations  (l4)  were  computed  for  boundary  conditions 
of  Q  ranging  over  three  orders  of  magnitude.  The 
solutions  can  be  applied  to  any  arbitrary  roughness  by 
accepting  as  the  boundary  condition  the  values  of  Q 
and  R  at  the  corresponding  height  in  the  solutions.  The 
roughness  plays  a  significant  role  in  translating  the 
solutions  of  Q  and  R  into  terms  of  real  physical 
variables.  By  definition,  we  know  that  the  real  velocity 
is  zero  at  the  roughness  height  z0.  In  terms  of  the 
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nondimensional  variables  this  means  that  u2  +  v2  =  1. 
This  can  be  expressed  according  to  equation  (10)  as 
2  3s  2 ,  3t  2  _  This  equation  can  be  used  to 
Ro  Tx  +  TZ  ~  1* 

determine  the  surface  Rossby  number,  R  ,  for  a  given 
solution  at  any  roughness  zQ  by  the  equation 


V) 

O 


(Vzo)2 


or  the  quantity 


/ 3Q2  +  V/2 

\3X  3  X  J 

G/no  =  (zc/^)3 


(15) 
2 

(16) 


From  a  practical  standpoint  it  was  found  that  the 
evaluation  Of  the  derivatives  of  Q  and  R  at  the  lower 
boundary  could  not  be  adequately  determined  from  finite 
difference  formulas.  Consequently,  an  approximate 
analytical  solution  was  developed.  Q  and  P  are  nearly 
invariant  with  height  near  the  lower  boundary. 
Therefore,  we  can  express  Q  and  P  in  terms  of  their 
average  value  between  adjacent  grid  points  in  the  first 
meter  above  the  lower  boundary.  Furthermore,  we  know 
from  equation  (6)  that  y  -  (z  /L)x  near  the  boundary. 
Consequently,  the  first  equation  of  equations  (1^1) 
can  be  approximately  expressed  as 


3 2  R  _  i 

3x2  10 1 2 Pzn  x 

°  (17) 

where  the  bar  designates  average  values.  If  we 
integrate  equation  (17),  we  obtain  (let  _  A) 

1012Pzo 

3 R  _  3Flj  +  A  In  x 

3X  ~  3 X  (18) 

where  the  subscript  1  refers  to  the  lower  boundary  at 
which  x  =  1.  Integrating  equation  (18)  we  obtain 

dR, 

R  =  Rj  +  dx  ( x  -  1 )  +  A  ( x  In  x  **  x  +  1 ) 


(19) 


If  R  and  x  correspond  to  the  first  grid  point  above 
the  boundary ,  then  the  derivative  at  the  boundary  can 
be  evaluated  by 


dR i  =  R?  -  Ri 
dx  x2  -  1 


A 


‘x2  In  x2 
(x2  -  1) 


1 


(20) 


where  P  and  Q  are  the  average  values  of  Q  and  P  in  the 
layer  from  x  =  1  to  x2* 

Using  equation  (19)  we  can  express  the  second  of 

equations  (14)  as  (let  _L _  =  3) 

10  12Pzc 


iii 

ax2 


B 


Si  +  iSi  (1  -  i)+  A  iln  x  +  (7 
X  3x  X  l  x 


(21) 


Integrating  equation  (21)  we  obtain 

~ )  In  x  +  ;  {x  In  x 
ax  ax  L  \  ax  / 

+(H‘  -  A)(x  - i:] 

Integrating  equation  (22)  we  obtain 


(x  - 

(22) 


1)1 


-1} 


(23) 


Usi 

the 


ng  the  values  of  -1  and  x  at  the  second  grid  point, 
first  derivative  can  be  evaluated  by 


^X2.-ln_x2_ 


(24) 


The  second  tern  on  the  right  of  equation  (24)  is 
generally  nearly  trivial  in  comparison  with  the  first 
term.  The  numerical  values  of  the  derivatives  obtained 
from  equations  (20)  ana  (24)  can  be  used  to  evaluate 
the  derivatives  at  any  arbitrary  height  x  from  equations 
(18)  and  (22).  Therefore,  the  Rossby  number  and  G/f 
can  be  determined  for  any  arbitrary  roughness  length 
from  equations  (15)  and  (16).  Table  I  presents  the 
values  .of  the  Rossby  number  and  G/f lQ  for  various 
roughness  lengths  as  determined  from  the  various  solu¬ 
tions.  It  is  to  be  noted  that  for  each  solution  G/f l 
decreases  with  increasing  roughness, 

A  graphical  representation  of  u#/G  as  a  function 
of  G/f l  and  z0/L0  is  presented  in  Figure  2.  The 
quantity  u*/3  =  C  is  analogous  to  a  drag  coefficient. 
From  Figure  2  the  magnitude  of  u#/G  can  be  determined 
for  any  combination  of  G/f£  and  zQ/l  .  Thus,  the 
surface  shearing  stress  canoe  determined  when  the 
gecstrophic  wind  is  known;  that  is. 


o  C2  G2 

between  the  surface  -wind  and  the 


(25) 


r.d  can  be  evaluated  by 


s  ft  —  ,  4-  ft  3  <*■ 

ax 


/  !Li  -  arctan  (R/~) 
'  ax  < 


where  the  derivatives  are  referred  to  the  roughness 
height . 


TABLE  I 

Values  of  the  Surface  Rossby  Number  and  G/f lQ  for  Various 
Roughness  Lengths  from  Three  Solutions 
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RELEVENT  SOLUTION 


Figure  3  represents  a  graphical  representation 
of  a  as  a  function  of  G/f£  and  zQ/l0. 

According  to  equations  (11)  the  velocity  defects 
(U,  V)  are  proportional  to  the  derivatives  of  the 
shearing  stress  components.  These  derivatives  were 
determined  from  the  approximate  analytical  solution 
near  the  surface  and  the  remaining  derivatives  at  grid 
points  were  approximated  from  second-order  Lagrangian 
polynomials  using  values  at  adjacent  grid  points. 

Once  the  derivatives  are  determined  and  the  surface 
Rossby  number  is  computed  from  equation  (15)  the  vertical 
distribution  of  the  velocity  defects  can  be  comp'.. ted. 

From  the  denned  variables  the  normalized  stresses, 
friction  velocity,  and  eddy  viscosity  distributions 
can  also  be  found.  These  parameters  were  tabulated 
for  each  of  the  solutions. 


II/.  BOUNDARY  LAVE  R  HOVc  L 


To  this  point  we  have  been  purely  formal  in  that 
we  have  stated  the  requisite  differential  equations, 
postulated  a  particular  functional  form  for  the  mixing 
length, and  then  proceeded  to  find  a  number  of  solutions 
of  the  equations  using  the  assumed  mixing-length 
relationship.  The  next  important  question  is  Which  of 
these  solutions  are  relevant  to  the  real  atmospheric 
planetary  boundary  layer?  In  attempting  to  answer  this 
question  we  shall  invoke  the  important  notion  of 
similarity.  Briefly,  we  assume  that  cSie  Reynolds 
number  of  the  planetary  boundary  layer  is  so  large 
that  there  exists  a  universal  equilibrium  structure  for 
the  turbulent  eddies,  and  that  the  scale  of  the  turbulent 
eddies  as  measured  by  the  mixing  length  is  universally 
related  to  the  gross  scale  (depth)  of  the  boundary  layer 
for  any  given  Rossby  number,  It  remains  to  determine 
what  scale  parameters  are  to  be  used  to  measure  the 
eddy  sizes  and  boundary  layer  depth.  It  is  immediately 
apparent  from  equation  (10)  that  either  iQ  or  L  is 
adequate  to  characterize  the  eddy  structure;  however, 
the  selection  of  a  characteristic  scale  length  for  the 
boundary  layer  is  somewhat  more  difficult.  The 
following  lists  some  of  the  possible  scale  heights:  • 

-  G/f 

H2  =  U#/f 

H3  =  h  (gradient  w ind  level) 

=  (G/f)  C2 
il5  -  (G/f)  C2  C03 c* 

H6  ■  (G/f)  C2  sino 

A  re  iabibhshipi  slich  as  h j  was  Used  ih  the  first  paper 
by  Ulaokadar4,  but  his  second  report5  used  a  relation 
that  can  be  reduced  to  lls,  11 3  in  a  simple  and  direct 
scale  length  analogous  to  the  boundary  layer  depth  in 
fluid  mechanics.  Hi,  is  the  ratio  of  surface  stress  to 
pressure  gradients,  i.e.,  t  /(3p/3y).  H5  is  an 
integral  scale  which  can  beshown  to  be  equivalent  to 
Hs  =  /“  (v/G)dz,  that  is,  the  integral  of  the  normalized 
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cros s— i s c c ari c  f  1  ov,' ,  .his  scale  heirht  h a s  the 
interestin'*  property  that  5K ,  is  the  total  mass 
transfer  across  the  isobars.  H - ,  another  integral 
scale,  is  equivalent  to  (u/'")dn,  the  interral  of 
the  iscbaric  flov;  and  o  1;: *  is  the  total  mass  transport 
alcnr  the  isobars.  There  is  no  a_  priori  reason  ‘’or 
selecting  one  scaii’nr:  length  .over  another.  Hence  a 
comparison  of  tv-c-  basic  parameters  that  result  from 
the  model  ■.-.•as  made  '"or  each  cf  the  scaling  lengths. 

These  are  the  relationships  between  ti'°  an~le  (a) 

-'rj  ^ on  ^^CccrcC'^i  '*  6*^  ~  j  r,oi^r, ®c  ti *'F  ^ 

the  surface  Rossby  number. 

Hj  tends  to  result  in  higher  than  observed  dram 
coefficients  at  high  surface  P.cssby  numbers,  although 
a's  are  suite  small  as  observed.  H2  and  H3  rive 
nearly  the  same  results  and  rive  more  reasonable  drag 
coefficients,  but  slightly  larger  angles.  Hu  and 
h 5  again  rive  comparable  results.  The  drag  coefficients 
are  more  in  agreement  with.  observations,  but  the  a '  s 
again  are  larger  than  observed.  H6  rives  a  drar 
coefficient  in  closer  agreement  ::ith  observations,  but 
has  the.  ooorest  agreement  for  a's. 


he  have  selected  H3  as  our  characteristic  ler.rth 
scale  because  cf  its  simplicity  and  because  preliminary 
results  show  its  adaptability  to  the  thermal  wind  case. 

According  to  our  notion  of  similarity,  the  nuar. 

( ii 3/ £  ),  representing  the  ratio  of  the  length  scale  0 
the  boundary  layer  and  the  length  scale  of  the  eddy 
structure,  is  a  universal  constant.  To  determine  its 
magnitude  we  resort  to  the  data  from  the  Leipzir  wind 
profile.  These  data  are  ri3  =  1071m  and  l0  =  32m  with 
the  result  that  K3/t_  =  33*  h'ith  this  constant  we 

are  able  to  select  tr.e  solution  which  is  o reposed  as  the 
universal  solution  of  the  steady  state,  barotrcpic, 
neutral  planetary  boundary  layer.  Table  II  presents  this 
solution  in  terms  of  the  r.or.dinensional  variables. 


be  can  new  demonstrate  how  tr.is  solution  is  used  to 
determine  the  distribution  cf  wind  and  turbulence  in 
the  planetary  boundary  layer.  ,'e  can  also  illustrate 
some  salient  features  cf  tne  model.  To  use  the  model, 
it  is  r.° ~ess ary  to  specify  the  three  parameters  which 
determine  the  surface  Rcsscy  number— namely ,  the 
reostrcphic  wind,  the  latitude,  and  the  surface  rourhness. 
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The  geostrophic  wind  can  be  readily  evaluated  from 
synoptic  (or  prognostic)  weather  maps,  while  several 
sources  6-7  can  be  consulted  to  evaluate  the  roughness 
at  the  site. 

The  relationship  of  Hf/G  and  £0f/G  to  the  surface 
Rossby  number  as  derived  from  the  model  is  illustrated  in 
Figure  4.  The  magnitude  of  l0  and  H  ean  be  determined 
from  Figure  4;  consequently,  the  nondimensional  height 
(z ,/v.)  of  Table  II  can  be  converted  to  true  height  z. 

Details  of  the  computation  of  this  and  other  dimensional 
parameters  for  the  boundary  layer  are  given  in  Appendix  A. 
From  Figure  4  It  can  be  seen  that  for  a  fixed  (G/f),  an 
increase  in  roughness  Increases  the  depth  of  the  boundary 
layer  (H)  and  the  limiting  mixing  length  (£u).  In  a 
like  manner,  for  a  fixed  roughness,  an  increase  in  the 
geostrophic  wind  increases  (H)  and  (tQ).  Decreasing  the 
latitude  tends  to  increase  these  parameters.  It  should 
be  noted  that  the  boundary  layer  depth  is  rather  insensitive 
to  roughness. 

We  can  also  determine  the  relationship  of  the  angle 
of  the  surface  wind  and  the  geostrophic  drag  coefficient 
to  the  surface  Rossby  number.  The  model  relationships 
are  illustrated  by  the  curves  of  Figure  5.  In  addition, 
all  known  observations  (Table  III)  have  been  plotted  to 
permit  comparison  with  the  results  of  the  model.  Small 
angles  of  the  surface  wind  are  associated  with  strong 
geostrophic  winds,  low  latitude,  and  small  roughness.  The 
trend  of  the  geostrophic  drag  coefficient  Is  similar  to 
that  of  the  surface  wind.  Unlike  fully  turbulent  flow 
in  ducts,  the  drag  coefficient  of  the  planetary  boundary 
layer  is  not  independent  of  the  flow  velocity,  thus, 
the  surface  shearing  stress  as  evaluated  from  equation  (19) 
is  not  proportional  to  the  square  of  the  geostrophic 
wind,  but  to  a  power  less  than  two.  Comparing  Figure  4 
with  Figure  5,  It  is  apparent  that  the  surface  wind  angle 
and  the  geostrophic  drag  coefficient  are  less  sensitive 
to  changes  in  G/f  than  are  H  and  l0.  With  the  exception 
of  the  Leipzig  wind  profile,  the  available  data  do  not 
satisfy  the  conditions  of  the  model  (barotropic,  neutral, 
steady  state)  or  the  method  of  computation  used  to  analyze 
the  data  is  of  doubtful  validity."  There  is  qualitative 
agreement  between  the  model  and  the  observations,  but  the 
scatter  is  too  great  and  the  observations  are  too  few  to 
permit  a  formal  statistical  analysis  of  the  results. 
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Graphic  presentations  of  salient  features  of  the 
model  are  presented  in  Figures  6  through  9-  Figure  6  is 
a  hodograph  of  the  wind  vector  in  arbitrary  units.  To 
use  the  graph  a  line  is  drawn,  as  illustrated  in  the 
figure,  from  the  origin  to  the  point  on  the  hodograph 
corresponding  to  *-he  value  of  z/£Q.  This,  line  represents 
the  geostrophic  wind  vector.  The  velocity  at  any 
arbitrary  height  z  is  determined  by  drawing  a  line  from 
the  point  on  the  hodograph  corresponding  to  z  /£  to  the 
desired  z/ Z  In- a  similar  manner,  the  hodograph  of 
the  shearing  stress  vector  is  shov/n  in  Figure  7.  To 
determine  the  magnitude  of  the  surface  stress,  we  obtain 
the  geostrophic  drag  coefficient  from  Figure  g b  and  then 
use  equation  (19). 

Figure  8  shows  the  vertical  distribution  of  the 
relative  eddy  viscosity  k/£.QG.  From  orevious  discussion 
we  concluded  that  for  a  fixed  roughness  and  latitude, 
l  increases  with  increasing  geostrophic  wind.  Therefore, 
the  eddy  viscosity  is  proportional  to  the  geostrophic 
wind  to  a  power  greater  than  one.  A  simple  computation 
shows  that  the  maximum  of  eddy  viscosity  occurs  at  a 
height  approximately  one-fifth  of  the  gradient  level. 
Figure  9  is  a  graphic  representation  of  the  vertical 
distribution  of  the  relative  viscous  dissipation  (e£0/G3), 


FIG.  6  WIND  VELOCITY  HODOGRAPH 

Best  Available  Copy 
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SURFACE  TOiA0*I.OXI 


V.  V1SCUSS10H 


The  planetary  boundary  layer  model  which  has 
been  presented  contains  the  features  sought  in  the 
objective.  In  the  ^ee  atmosphere,  the  wind  approaches 
the  geostrophic  wind,  while  in  the  surface  boundary 
layer  the  vertical  distribution  of  wind  is  logarithmic. 
However,  it  would  be  misleading  to  imply  from  this  that 
the  idealized  planetary  boundary  layer  model  could  be 
as  successful  as  these  two  otherwise  independent  concepts. 
The  success  of  the  geostrophic  wind  eouation  lies  in  the 
fact  that  the  wind  appears  to  respond  quickly  to  pressure 
gradient  changes,  with  the- result  that  the  wind  is 
quas.i-geostrophic  in  the  free  atmosphere  even  when  condi- 
tions  are  not  ideal.  In  like  manner,  the  wind  profile 
responds  quickly  to  stress  and  roughness  changes  in  the 
surface  boundary  layer.  However,  Batchelor  (Reference  8, 
o.  2R6)  has  emphasized  that  the  stress  at  the  ground  is 
likely  to  respond  to  pressure  gradient  changes  slowly, 
with  the  result  that  a  true  steady-state  planetary 
boundary  layer  would  be  a  rare  phenomenon  in  the  real 
atmosphere.  In  reply,  Lettau  (Reference  8,  p.  257)  has 
stated  that  in  the  lower  portion  of  the  boundary  layer 
the  profile  will  respond  relatively  quickly  to  pressure 
gradient  changes,  but  more  prolonged  inertial  oscillations 
will  complicate  the  upper  portion  of  the  boundary  layer. 

In  the  course  of  solving  equations  (13),  the  authors 
met  a  problem  of  computational  instability  in  the  upper 
portion  of  the  boundary  layer  solution.  This  Instability 
had  all  the  appearances  of  inertial  oscillations  of  the 
type  expected  by  Lettau.  However,  in  spite  of  these 
oscillations,  the.  lower  portion  of  the  profile  was 
stable  and  the  distribution  of  stress  was  in  consonance 
with  the  steady-state  solution.  This  result  is  rather 
encouraging;  it  implies  that  at  the  boundary  such 
features  as  the  stress  and  angle  of  the  wind  are  in 
balance  with  the  geostrophic  wind ,  even  if  the  upoer 
portion  of  the  boundary  layer  remains  in  oscillation. 

In  any  event,  the  assumption  of  steady-state  sets  a 
restriction  on  the  utility  of  the  boundary  layer  solutions. 
A  quantitative  evaluation  of  this  limitation  ’would 
require  that  the  "time  constant"  of  the  boundary  layer  be 
evaluated  through  non-steady  solutions. 
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The  restriction  of  the  model  to 
is  another  obvious  limitation  of  the 
state  neutral  boundary  layer  model  is 
applicable  to  clear  sky  conditions  cv 
it  should  be  often  applicable  tc  core 
or  over  i and  surfaces  with  evert -st  s 


a  ci  i  nb  a  t  i  c  c  o  n  d  i  t  i  -  n  s 
model.  A  5teafl.v- 
r.ot  like  lv  to  be 
er  land.  However , 
it  ions  ever  the  sea 
\*  ~  0  ~  , 


The  restriction  of  the  solutions  to  barotroni c 
conditions  war  not  essential,  but  it  did  greatly 
simplify  the  problem  in  terms  of  the  number  of  solutions 
required.  Furthermore,  it  would  seem  that  for  the  case 
of  constant  thermal  wind,  the  boundary  layer  model  in 
terms  of  the  deviation  Gf  the  wind  from  the  geostroohic 
wind:  would  remain  valid.  Consequently,  except  for  extreme 
thermal  wind,  it  is  likely  that  the  barcolinic  model  can 
be  approximated  by  inclusion  of  thermal  wind  in  the 
barotronic  model.  The  validity  of  this  postulate 
remains  a  matter  to  be  tested. 


It  is  not  too  difficult  to  visualize  the  limitation 
of  the  model  due  to  the  assumptions  of  steady-state 
oarotropic  and  adiabatic  conditions.  Aside  from  these, 
the  critical  issue  regarding  the  planetary  boundary  layer 
model  is  the  assumption  involved  in  expressing  the  stress 
in  terms. of  the  mean  velocity.  Here  it  was  necessary  to 
introduce  the  concepts  of  the  mixing-length  hypothesis. 

In  some  circles,  the  mixing-length  hypothesis  is  in 
disfavor  since  it  is  based  more  on  intuition  than 
fundamentals.  nevertheless,  the  mixing-length  approach 
has  beer,  successful  in  some  areas  of  turbulence  and  it 
remains  today  the  only  oractical  approach  for  prediction 
mode  Is . 
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VI.  COMPARISON  ulITH  OTHER  HOVELS 


Similar  models  based  on  the  mixing- length  hypothesis 
have  been  reported  by  Blackadar4"5  and  Lettau9.  Each 
model  is  based  on  a  slightly  different  mixing  length  and 
the  equa‘ ' ons  were  solved  by  different  numerical  methods. 
Lettau  used  a  mixing  length  derived  from  duct  flow  and 
Blackadar  used  an  empirical  mixing  length  similar  in 
form  to  the  one  used  in  this  study. 

A  comparison  of  the  drag  coefficients  and  surface 
wind  as  functions  of  surface  Rossby  numbers  from 
the  three  models  is  shown  in  Figures  10a  and  b.  It  can 
be  seen  that  the  three  models  give  similar  results. 
Comparing  these  results  with  the  observation  on  Figure  5 
indicates  that  one  model  cannot  be  singled  out  as  being 
in  better  agreement  with  observations  than  the  other. 


FIG  10  A  COMPARISON  OF  THE  SURFACE  WIND  ANGLES  AND  THE  DRAG 
COEFFICIENTS  PREDICTED  BY  THE  THREE  -  MODELS. 
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VII. 


CONCLUSIONS 


We  believe  that  the  development  and  solution  of  these 
planetary  boundary  layer  models  represents  a  significant 
step  forward,  since  it  furnishes  a  connecting  link 
between  mic’-'orneteorology  and  macrometeorology.  We  do 
not  contend  that  the  model  presented  is  the  final  answer 
as  it  is  based  on  several  hypotheses  which  require 
experimental  verification.  The  model  is  restricted  to 
a  class  of  atmospheric  conditions.  However,  it  represents 
a  foundation  which  can  be  refined  as  more  experimental 
data  becume  available,  and  extended  to  cover  a  broader 
class  of  meteorological  conditions. 

The  present  model  has  strategic  and  design  applica¬ 
tions  which  require  a  knowledge  of  the  distribution  of 
wind  or  turbulence  characteristics  for  differing  situations. 
For  example,  it  is  believed  that  the  model  can  be  used 
to  evaluate  the  effect  of  roughness  on  the  diffusion  of 
CUR  agents.  Such  a  study  is  now  under  way.  In  like 
manner,  the  model  can  be.  extended  to  a  broader  .class  of 
micrometeorological  problems  such  as  determining  climatic 
and  synoptic  estimates  of  evapotranspiration  from  large- 
scale  watersheds,  and  turbulent  transfer  of  heat,  water 
vapor,  and  carbon  dioxide  from  the  earth's  surface. 

From  a  long  range  standpoint,  one  of  the  great 
challenges  facing  meteorology  is  the  development  of  a 
truly  diabatic  model  for  the  prediction  of  the  hemispheric 
circulation.  An  important  feature  of  such  a  model  would 
be  the  energy  dissipation  associated  with  the  planetary, 
boundary  layer.  The  model  furnishes  estimates  of  the 
dissipation  and  surface  shearing  stress  derived  from  • 
knowledge  of  only  the  geostrophic  wind,  latitude  and 
roughness. 

It  is  concluded,  therefore,  that  the  planetary  boundary 
layer  model  presented  has  useful  applications  and  is  worthy 
of  extension  to  cover  a  broader  class  of  situations. 
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APPENDIX  A 


Computation  oj  VinuinAional  Van.amzte.Ki>  oj  the,  Boundary  Lai/e.A 


In  order  to  compute  dimensional  parameters  of  the 
boundary  layer,  the  geostrophic  wind,  latitude  (Coriolis 
parameters)  and  surface  roughness  for  the’  particular 
situation  must  be  determined.  Then  the  surface  Rossby 
number  and  G/f  can  be  computed.  From  Figure  ^  the  limiting 
mixing  length  (iQ)  and  the  boundary  layer  depth  can  be 
determined.  The  real  height  z  can  then  be  obtained  from 
the  z / lQ  column  of  Table  II.  The  real  wind  components  u 
and  v  for  each  level  can  be  computed  in  two  ways  depending 
on  the  orientation  desired.  In  the  original  orientation 
with  the  surface  winds  along  the  x  axis,  the  following 
equations  are  applicable: 

u  =  [cosa0  -  r/rQ  cos  a  ]G  and  v  =  [sinaQ-  r/rQsin  ct]G 

If  one  desires  to  orient  the  geostrophic  wind  along  the 
x  axis,  the  equations  are 


u  =  [T  -  r/rQ  cos  (a-a0)]G 


v  =  [r/rQ  sin  (a-a0)]G 
where  r0  is  the  length  of  /gx-u\  2  +  /Gy-v\ 

\  G  /  z0  V5"”  /  z0 
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at  the  x0/£0  height  and  r  is  this  quantity  at  specific 
levels;  a  Is  the  angle  between  the  surface  and  geostrophic 
wind;  a  =  arctan  (Gy-v)/G/(Gy-v)/G.  To  simplify  the 
computation,  value-  bf  r  and  a  as  functions  of  z/l^ 
are  given  in  Table  IV,  Figure  11  shows  £„  as  a  function 
of  surface  Rossby  number;  a  can  be  found  from  Figure  10a. 


The  shearing  stress  components  are  obtained  in  a 
similar  manner.  First,  the  drag  coefficient  U#/G  must  be 
determined  from  Figure  10b;  tq  is  then  obtained  from  the 
equation  tq  «=  p  2  G2  ;  Tx°and  ry  can  be  found  from 

the  following  equations  depending  on  the  desired  orientation. 
For  the  surface  wind  oriented  along  the  x  axis,  the 
equations  are 

T  3  (t/tQCOS'l')T 
X  w  o 

Ty  =  (t/t0sin?)T0 
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Table  IV 

r,  t,  a  and  y  as  functions  of  z/lQ 


z/lo 

r 

a 

txlO3 

Y 

53.1 

0 

"A 

0 

Ct 

50.0 

.0023 

382.2 

.0009 

324.2 

45.0 

.0067 

309.4 

.0035 

261.9 

40.0 

.0150 

258.5 

.0093 

211.0 

35.0 

.0207 

237.2 

.0152 

187.5 

30.0 

.0466 

183.7 

.0449 

134.3 

25.0 

.0721 

154.1 

.0802 

105.2 

20.0 

.1054 

129.1 

.1332 

80.8 

15.0 

.1475 

108.8 

.2091 

58.4 

10.0 

.2013 

87.2 

.3144 

38.2 

6.4 

.2485 

73.4 

.4117 

24.7 

5.0 

.2699 

68.0 

.  4549 

20.0 

2.5 

.3195 

57.4 

.5418 

10.5 

1.0 

.3717 

48.6 

.5999 

04.8 

0.5 

.  4085 

43.6 

.6206 

2.5 

0.25 

.  4442 

39.7 

.6320 

1.5 

0.10 

.  4908 

35.4 

.6383 

0.3 

.01 

.6440 

26.6 

.6413 

0 

.001 

.7910 

21.2 

.6426 

0 

.0001 

.9300 

18.0 

.6428 

0 

.00001 

1.069 

15.5 

.6428 

b 

.000001 

1.209 

13.7 

.6428 

0 

01  X 


For  the  geostropnic  wind  oriented  along  the  x  axis 


Tx  =  t/toCOSUo-^TQ 
ty  =  t/t0sin(a0-Y)i0 

where  t  =  [(tx/pG2)2  +  { ty/pG2 ) 2 ] 2  ;  t  as  a  function  of 
z /l  is  given  in  Table  IV  and  tQ  is  given  in  Figure  11  as 
a  function  of  surface  Rossby  number;  ¥  =  arctan 
( ty/pG2  )z/(  -rx/pG2 )  z  is  also  given  in  Table  IV  as  a  function 
of  z/lQ  ;  =  0  and  a0  is  the  angle  of  the  surface  wind 

with  the  geostropnic. 

The  friction  velocity  U#  is  a  scalar  and  can  be 
found  by  multiplying  the  column  U#/G  by  a  constant.  The 
simplest  method  of  determining  the  constant  is  to  compute 
the  ratio  of  the  U*/G  given  in  Figure  10  for  the  desired 
surface  Rossby  number  to  the  (U*/G)  in  Table  II  for  the 
required  z Q/ 10  height.  This  ratio  must  be  multiplied  by 
G,  the  geostrophlc  wind  speed. 

Thus  U#  =  A  (iy 

\G  Az) 

wnero  A  =  (U*/G)„  /(U#/G)  ..  G 

no  o'  o 

To  determine  the  distribution  of  K  with  height,  the 
value  of  K/£qG,  found  in  the  last  column  of  Table  II  must 
also  be  multiplied  by  a  constant.  Since  K  =  U#£,  the 
above  constant  (A)  times  t0  is  used. 

Thus  K  =  ( K/ laG )  hi 

^  z  o 

A  sample  computation  for  all  parameters  of  the  Leipzig  wind 
profile  follows 

geostrophlc  wind  (G)  17.51  m/sec 

surface  roughness  (zQ)  .20  m 

Coriolis  parameter  (f)  1.14  x  10“1<sec“1 

air  density  (o)  1.15  x  10“3 

R0  *  G/f  z0  =  7.75  x  10 5;  log  R0  »  5.89 
G/f  =  1.538  x  105 
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To  compute  real  height  and  depth  of  the  boundary 
layer  from  Figure  (Hf/G)  x  103  =  7.0  ;  (£Qf/G)  x  104=  2.1 

H  =  1075;  l0  =  32.3  ;  .z Q/l0  =  6.19  x  10~3 

To  compute  u,  v,  tv  and  t 

x  y 

from  Figure  10  aQ  =  25.5  and  U*/G  =  .039  =  C 
from  Equation  25  =  pC2G2  =  5.36  dynes/cm2 

from  Figure  11  r^  =  .633  ;  tQ  =  0.64253  x  10“3 
To  compute  U* 

U*  =  Mtt2-)  where  A  =  (U»/G)R  /(U*/G)  ..  G 

v  /  Z  nO  Z0/ £q 

from  Table  II,  U#G  at  z Q/lQ  of  6.39  x  10"3  =  .0253 

A  =  .039  x  17.51 

- 70253 —  ~  27 

U#  =  27  X  .0253  =  .683  m/sec 

To  compute  K 

K  =  (K/£QG)z 

=  .00009  x  27  x  32.3  =  .078  mZ/sec 

Using  these  values.  Table  II  and  Table  IV  give  the  results 
shown  in  Table  V.  Figures  12  and  13  give  a  comparison  of 
parameters  derived  from  the  model  with  those  computed  by 
Lettau  3» 10  from  the  Leipzig  profile. 
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Leipzig  Wind  Profile  Parameters  Computed  from  the  Model 
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Computed  for  a  westerly  geostrophic  wind 


PARAMETER  DERIVED  FROM  OBSERVATIONS 
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The  objective  of  this  study  is  to  develop  a  theoretical 
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this  relationship  in  the  equation  of  motion  led  to  a  set  of  second 
order,  nonlinear  differential  equations,  which  were  solved  on  a 
digital  computer.  Universal  profiles  of  the  wind,  stress,  and  eddy 
viscosity  were  fixed  by  invoking  the  important  notion  of  similarity, 
that  is,  it  is  assumed  the  scale  of  turbulence  is  uniquely  related  to 
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constant  is  evaluated  from  experimental  data.  Possible  applications 
of  the  model  to  practical  problems  are  outlined. 


DD 


FORM 

1  JAN  94 


1473 


UNCLASSIFIED 


Security  Classification 


LINK  A 


UNK  B 


LINK  C 


UNCLASSIFIED 

Security  Classification 


key  wo  bos 


& 

Micrometeorology 

Numerical  Solutions 

Planetary  Boundary  Layer 

Wind  and  Turbulence  Distribution 

Dissipation 

Eddy  Viscosity 

Mixing  Length 


I 


INSTRUCTIONS 


1.  ORIGINATING  ACTIVITY:  Enter  the  name  and  addreaa 
of  the  contractor,  subcontractor,  grantee,  Department  of  De¬ 
fense  activity  or  other  organization  (corporate  oothor)  issuing 
the  report. 

2a.  REPORT  SECURITY  CLASSIFICATION:  Enter  the  over¬ 
all  security  classification  of  the  report.  Indicate  whether 
“Restricted  Data”  is  included.  Marking  ia  to  be  in  accord¬ 
ance  with  appropriate  security  regulations. 

2b.  GROUP:  Automatic  downgrading  is  specified  in  DoD  Di¬ 
rective  5200.10  and  Armed  Forces  Industrial  Manual.  Enter 
the  group  number.  Also,  when  applicable,  show  that  optional 
markings  have  been  used  for  Group  3  and  Group  4  as  author¬ 
ized. 

3.  REPORT  TITLE:  Enter  the  complete  report  title  In  all 
capital  letter  a.  Titles  in  ail  caaea  should  be  unclassified, 
if  a  meaningful  title  cannot  be  selected  without  classifica¬ 
tion.  show  title  classification  in  all  capitals  in  parenthesis 
immediately  following  the  title. 

4.  DESCRIPTIVE  NOTES:  If  appropriate,  enter  the  type  of 
report,  e.g.,  interim,  progress,  aummsry,  annual,  or  final. 

Give  ihe  inclusive  dates  when  a  specific  reporting  period  la 
covered. 

5.  AUTHOR(S):  Enter  the  narnefs)  of  authorfm)  aa  shown  on 
or  in  the  report.  Enter  last  name,  first  name,  middle  initial. 

If  military,  show  rank  and  branch  of  service.  The  name  of 
the  principal  author  is  an  absolute  minimum  requirement. 

6.  REiPORT  DATE:  Enter  the  date  of  the  report  as  day, 
month,  year,  or  month,  year.  If  more  than  one  date  appears 
on  the  report,  use  date  of  publication. 

:  7a  TOTAL  NUMBER  OF  PAGES:  The  total  page  count 
should  follow  normal  pagination  procedures,  he.,  enter  the 
number  of  pages  containing  information 

7b.  NUMBER  OF  REFERENCES  Enter  the  total  number  of 
references  cited  in  the  report. 

8  a.  CONTRACT  OR  GRANT  NUMBER:  If  spproprlate,  enter 
the  applicable  number  of  the  contract  or  grant  under  which 
the  report  waa  written. 

36.  8c.  8>  Bd.  PROJECT  NUMBER:  Enter  the  appropriate 
military  department  identification,  such  aa  project  number, 
subproject  number,  system  numbers,  task  number,  etc. 

9a.  ORIGINATOR'S  REPORT  NUMBER(S):  Enter  the  offl. 
cial  report  number  by  which  the  document  will  be  identified 
and  cortrolted  by  the  originating  activity.  This  number  must 
be  unique  to  thli  report. 

96.  OTHER  REPORT  NUMBERfS):  If  the  report  has  been 
assigned  sny  other  report  numbers  (eithar  by  the  originator 
or  by  the  sponsor),  also  enter  this  numberfa). 


10.  AVAILABILITY/LIMITATION  NOTICES:  Enter  any  lim¬ 
itations  on  further  dissemination  of  the  report,  other  than  those 
imposed  by  security  classification,  using  standard  statements 
such  is: 

(1)  "Qualified  requesters  may  obtain  copies  of  this 
report  from  DDC.” 

(2)  "Foreign  announcement  and  dissemination  of  this 
report  by  DDC  is  not  authorized." 

(3)  "U.  S.  Government  agencies  may  obtain  copies  of 
this  report  directly  from  DDC.  Other  qualified  DDC 
users  shall  request  through 

ft 

(4)  "U.  S-  military  agencies  may  obtain  copies  of  this 
report  directly  from  DDC.  Other  qualified  uaera 
shall  raqueat  through 

(3)  "All  distribution  of  this  report  is  controlled.  Qual¬ 
ified  DDC  users  shall  request  through 

*t 

If  tha  report  has  been  furnished  to  the  Office  of  Technical 
Services,  Department  of  Commerce,  for  sale  to  the  public,  indi¬ 
cate  this  fact  and  enter  the  price,  if  known. 

It  SUPPLEMENTARY  NOTES:  Use  for  additional  explana¬ 
tory  note* 

12.  SPONSORING  MILITARY  ACTIVITY:  Enter  the  name  of 
the  dopsrtmontal  project  office  or  laboratory  sponsoring  (pay- 
ini  lor)  iho  research  and  development.  Include  address- 

13.  ABSTRACT:  Enter  an  abstract  giving  a  brief  and  factual 
summary  of  tha  document  indicative  of  the  report,  even  though 
it  may  alao  appear  elaewhere  in  the  body  of  the  technical  re¬ 
port,  If  additional  apace  ia  required,  a  continuation  aheet 
aball  be  attached, 

It  is  highly  desirable  that  the  abstract  of  classified  re¬ 
ports  be  unclassified.  Each  paragraph  of  the  abetract  ahall 
end  with  an  Indication  of  the  military  aecurity  Classification 
of  tha  information  in  the  peragraph,  represented  as  fTS),  (S), 
(C),  or  (V). 

There  la  no  limitation  on  the  length  of  the  abetract.  How¬ 
ever,  Ihe  euggeated  length  is  from  150  to  225  words. 

14.  REV  WORDS:  Kay  words  are  technically  mesoingful  terms 
or  abort  phrases  that  characterize  a  report  and  may  be  used  aa 
Index  entries  for  cataloging  the  report.  Key  words  must  be 
attested  aa  that  no  security  classification  is  required. '  Iden- 
fiars,  such  aa  equipment  model  designation,  trade  name,  mili¬ 
tary  project  coda  name,  geographic  location,  may  be  used  as 
key  word!  but  wilt  be  followed  by  an  indication  of  technical 
context,  The  aaaignment  of  links,  rules,  and  weights  is 
optional. 


y,;  v.i  i 


4  Ik  ‘  J  ii  4  <ia»£  I  'l/ 


UNCLASSIFIED' 

Security  Classification 


